Gene-based burden scores identify rare variant associations for 28 blood biomarkers.

BACKGROUND: A relevant part of the genetic architecture of complex traits is still unknown; despite the discovery of many disease-associated common variants. Polygenic risk score (PRS) models are based on the evaluation of the additive effects attributable to common variants and have been successfully implemented to assess the genetic susceptibility for many phenotypes. In contrast, burden tests are often used to identify an enrichment of rare deleterious variants in specific genes. Both kinds of genetic contributions are typically analyzed independently. Many studies suggest that complex phenotypes are influenced by both low effect common variants and high effect rare deleterious variants. The aim of this paper is to integrate the effect of both common and rare functional variants for a more comprehensive genetic risk modeling. METHODS: We developed a framework combining gene-based scores based on the enrichment of rare functionally relevant variants with genome-wide PRS based on common variants for association analysis and prediction models. We applied our framework on UK Biobank dataset with genotyping and exome data and considered 28 blood biomarkers levels as target phenotypes. For each biomarker, an association analysis was performed on full cohort using gene-based scores (GBS). The cohort was then split into 3 subsets for PRS construction and feature selection, predictive model training, and independent evaluation, respectively. Prediction models were generated including either PRS, GBS or both (combined). RESULTS: Association analyses of the cohort were able to detect significant genes that were previously known to be associated with different biomarkers. Interestingly, the analyses also revealed heterogeneous effect sizes and directionality highlighting the complexity of the blood biomarkers regulation. However, the combined models for many biomarkers show little or no improvement in prediction accuracy compared to the PRS models. CONCLUSION: This study shows that rare variants play an important role in the genetic architecture of complex multifactorial traits such as blood biomarkers. However, while rare deleterious variants play a strong role at an individual level, our results indicate that classical common variant based PRS might be more informative to predict the genetic susceptibility at the population level.

Analysis of 72,469 UK Biobank exomes links rare variants to male-pattern hair loss.

Male-pattern hair loss (MPHL) is common and highly heritable. While genome-wide association studies (GWAS) have generated insights into the contribution of common variants to MPHL etiology, the relevance of rare variants remains unclear. To determine the contribution of rare variants to MPHL etiology, we perform gene-based and single-variant analyses in exome-sequencing data from 72,469 male UK Biobank participants. While our population-level risk prediction suggests that rare variants make only a minor contribution to general MPHL risk, our rare variant collapsing tests identified a total of five significant gene associations. These findings provide additional evidence for previously implicated genes (EDA2R, WNT10A) and highlight novel risk genes at and beyond GWAS loci (HEPH, CEPT1, EIF3F). Furthermore, MPHL-associated genes are enriched for genes considered causal for monogenic trichoses. Together, our findings broaden the MPHL-associated allelic spectrum and provide insights into MPHL pathobiology and a shared basis with monogenic hair loss disorders.

Identification of de novo variants in nonsyndromic cleft lip with/without cleft palate patients with low polygenic risk scores.

BACKGROUND: Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants. METHODS: To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization. CONCLUSION: In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.

Genome-wide identification of disease-causing copy number variations in 450 individuals with anorectal malformations.

Anorectal malformations (ARM) represent a spectrum of rare malformations originating from a perturbated development of the embryonic hindgut. Approximately 60% occur as a part of a defined genetic syndrome or within the spectrum of additional congenital anomalies. Rare copy number variations (CNVs) have been associated with both syndromic and non-syndromic forms. The present study represents the largest study to date to explore the contribution of CNVs to the expression of ARMs. SNP-array-based molecular karyotyping was applied in 450 individuals with ARM and 4392 healthy controls. CNVs were identified from raw intensity data using PennCNV. Overlapping CNVs between cases and controls were discarded. Remaining CNVs were filtered using a stringent filter algorithm of nine filter steps. Prioritized CNVs were confirmed using qPCR. Filtering prioritized and qPCR confirmed four microscopic chromosomal anomalies and nine submicroscopic CNVs comprising seven microdeletions (del2p13.2, del4p16.2, del7q31.33, del9p24.1, del16q12.1, del18q32, del22q11.21) and two microduplications (dup2p13.2, dup17q12) in 14 individuals (12 singletons and one affected sib-pair). Within these CNVs, based on their embryonic expression data and function, we suggest FOXK2, LPP, and SALL3 as putative candidate genes. Overall, our CNV analysis identified putative microscopic and submicroscopic chromosomal rearrangements in 3% of cases. Functional characterization and re-sequencing of suggested candidate genes is warranted.

Next-generation phenotyping contributing to the identification of a 4.7 kb deletion in KANSL1 causing Koolen-de Vries syndrome.

Next-generation phenotyping (NGP) is an application of advanced methods of computer vision on medical imaging data such as portrait photos of individuals with rare disorders. NGP on portraits results in gestalt scores that can be used for the selection of appropriate genetic tests, and for the interpretation of the molecular data. Here, we report on an exceptional case of a young girl that was presented at the age of 8 and 15 and enrolled in NGP diagnostics on the latter occasion. The girl had clinical features associated with Koolen-de Vries syndrome (KdVS) and a suggestive facial gestalt. However, chromosomal microarray (CMA), Sanger sequencing, multiplex ligation-dependent probe analysis (MLPA), and trio exome sequencing remained inconclusive. Based on the highly indicative gestalt score for KdVS, the decision was made to perform genome sequencing to also evaluate noncoding variants. This analysis revealed a 4.7 kb de novo deletion partially affecting intron 6 and exon 7 of the KANSL1 gene. This is the smallest reported structural variant to date for this phenotype. The case illustrates how NGP can be integrated into the iterative diagnostic process of test selection and interpretation of sequencing results.

Reconstruction of the origin of the first major SARS-CoV-2 outbreak in Germany.

The first major COVID-19 outbreak in Germany occurred in Heinsberg in February 2020 with 388 officially reported cases. Unexpectedly, the first outbreak happened in a small town with little to no travelers. We used phylogenetic analyses to investigate the origin and spread of the virus in this outbreak. We sequenced 90 (23%) SARS-CoV-2 genomes from the 388 reported cases including the samples from the first documented cases. Phylogenetic analyses of these sequences revealed mainly two circulating strains with 74 samples assigned to lineage B.3 and 6 samples assigned to lineage B.1. Lineage B.3 was introduced first and probably caused the initial spread. Using phylogenetic analysis tools, we were able to identify closely related strains in France and hypothesized the possible introduction from France.

GenRisk: a tool for comprehensive genetic risk modeling.

SUMMARY: The genetic architecture of complex traits can be influenced by both many common regulatory variants with small effect sizes and rare deleterious variants in coding regions with larger effect sizes. However, the two kinds of genetic contributions are typically analyzed independently. Here, we present GenRisk, a python package for the computation and the integration of gene scores based on the burden of rare deleterious variants and common-variants-based polygenic risk scores. The derived scores can be analyzed within GenRisk to perform association tests or to derive phenotype prediction models by testing multiple classification and regression approaches. GenRisk is compatible with VCF input file formats. AVAILABILITY AND IMPLEMENTATION: GenRisk is an open source publicly available python package that can be downloaded or installed from Github (https://github.com/AldisiRana/GenRisk). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

GestaltMatcher facilitates rare disease matching using facial phenotype descriptors.

Many monogenic disorders cause a characteristic facial morphology. Artificial intelligence can support physicians in recognizing these patterns by associating facial phenotypes with the underlying syndrome through training on thousands of patient photographs. However, this 'supervised' approach means that diagnoses are only possible if the disorder was part of the training set. To improve recognition of ultra-rare disorders, we developed GestaltMatcher, an encoder for portraits that is based on a deep convolutional neural network. Photographs of 17,560 patients with 1,115 rare disorders were used to define a Clinical Face Phenotype Space, in which distances between cases define syndromic similarity. Here we show that patients can be matched to others with the same molecular diagnosis even when the disorder was not included in the training set. Together with mutation data, GestaltMatcher could not only accelerate the clinical diagnosis of patients with ultra-rare disorders and facial dysmorphism but also enable the delineation of new phenotypes.

Heterogeneity and excitability of BRAFV600E-induced tumors is determined by Akt/mTOR-signaling state and Trp53-loss.

BACKGROUND: Developmental brain tumors harboring BRAFV600E somatic mutation are diverse. Here, we describe molecular factors that determine BRAFV600E-induced tumor biology and function. METHODS: Intraventricular in utero electroporation in combination with the piggyBac transposon system was utilized to generate developmental brain neoplasms, which were comprehensively analyzed with regard to growth using near-infrared in-vivo imaging, transcript signatures by RNA sequencing, and neuronal activity by multielectrode arrays. RESULTS: BRAF  V600E expression in murine neural progenitors elicits benign neoplasms composed of enlarged dysmorphic neurons and neoplastic astroglia recapitulating ganglioglioma (GG) only in concert with active Akt/mTOR-signaling. Purely glial tumors resembling aspects of polymorphous low-grade neuroepithelial tumors of the young (PLNTYs) emerge from BRAFV600E alone. Additional somatic Trp53-loss is sufficient to generate anaplastic GGs (aGGs) with glioneuronal clonality. Functionally, only BRAFV600E/pAkt tumors intrinsically generate substantial neuronal activity and show enhanced relay to adjacent tissue conferring high epilepsy propensity. In contrast, PLNTY- and aGG models lack significant spike activity, which appears in line with the glial differentiation of the former and a dysfunctional tissue structure combined with reduced neuronal transcript signatures in the latter. CONCLUSION: mTOR-signaling and Trp53-loss critically determine the biological diversity and electrical activity of BRAFV600E-induced tumors.

Iron Deficiency Caused by Intestinal Iron Loss-Novel Candidate Genes for Severe Anemia.

The adult human body contains about 4 g of iron. About 1-2 mg of iron is absorbed every day, and in healthy individuals, the same amount is excreted. We describe a patient who presents with severe iron deficiency anemia with hemoglobin levels below 6 g/dL and ferritin levels below 30 ng/mL. Although red blood cell concentrates and intravenous iron have been substituted every month for years, body iron stores remain depleted. Diagnostics have included several esophago-gastro-duodenoscopies, colonoscopies, MRI of the liver, repetitive bone marrow biopsies, psychological analysis, application of radioactive iron to determine intact erythropoiesis, and measurement of iron excretion in urine and feces. Typically, gastrointestinal bleeding is a major cause of iron loss. Surprisingly, intestinal iron excretion in stool in the patient was repetitively increased, without gastrointestinal bleeding. Furthermore, whole exome sequencing was performed in the patient and additional family members to identify potential causative genetic variants that may cause intestinal iron loss. Under different inheritance models, several rare mutations were identified, two of which (in CISD1 and KRI1) are likely to be functionally relevant. Intestinal iron loss in the current form has not yet been described and is, with high probability, the cause of the severe iron deficiency anemia in this patient.

Variant profiling of colorectal adenomas from three patients of two families with MSH3-related adenomatous polyposis.

The spectrum of somatic genetic variation in colorectal adenomas caused by biallelic pathogenic germline variants in the MSH3 gene, was comprehensively analysed to characterise mutational signatures and identify potential driver genes and pathways of MSH3-related tumourigenesis. Three patients from two families with MSH3-associated polyposis were included. Whole exome sequencing of nine adenomas and matched normal tissue was performed. The amount of somatic variants in the MSH3-deficient adenomas and the pattern of single nucleotide variants (SNVs) was similar to sporadic adenomas, whereas the fraction of small insertions/deletions (indels) (21-42% of all small variants) was significantly higher. Interestingly, pathogenic somatic APC variants were found in all but one adenoma. The vast majority (12/13) of these were di-, tetra-, or penta-base pair (bp) deletions. The fraction of APC indels was significantly higher than that reported in patients with familial adenomatous polyposis (FAP) (p < 0.01) or in sporadic adenomas (p < 0.0001). In MSH3-deficient adenomas, the occurrence of APC indels in a repetitive sequence context was significantly higher than in FAP patients (p < 0.01). In addition, the MSH3-deficient adenomas harboured one to five (recurrent) somatic variants in 13 established or candidate driver genes for early colorectal carcinogenesis, including ACVR2A and ARID genes. Our data suggest that MSH3-related colorectal carcinogenesis seems to follow the classical APC-driven pathway. In line with the specific function of MSH3 in the mismatch repair (MMR) system, we identified a characteristic APC mutational pattern in MSH3-deficient adenomas, and confirmed further driver genes for colorectal tumourigenesis.

Genome-Wide Survey for Microdeletions or -Duplications in 155 Patients with Lower Urinary Tract Obstructions (LUTO).

Lower urinary tract obstruction (LUTO) is, in most cases, caused by anatomical blockage of the bladder outlet. The most common form are posterior urethral valves (PUVs), a male-limited phenotype. Here, we surveyed the genome of 155 LUTO patients to identify disease-causing CNVs. Raw intensity data were collected for CNVs detected in LUTO patients and 4.392 healthy controls using CNVPartition, QuantiSNP and PennCNV. Overlapping CNVs between patients and controls were discarded. Additional filtering implicated CNV frequency in the database of genomic variants, gene content and final visual inspection detecting 37 ultra-rare CNVs. After, prioritization qPCR analysis confirmed 3 microduplications, all detected in PUV patients. One microduplication (5q23.2) occurred de novo in the two remaining microduplications found on chromosome 1p36.21 and 10q23.31. Parental DNA was not available for segregation analysis. All three duplications comprised 11 coding genes: four human specific lncRNA and one microRNA. Three coding genes (FBLIM1, SLC16A12, SNCAIP) and the microRNA MIR107 have previously been shown to be expressed in the developing urinary tract of mouse embryos. We propose that duplications, rare or de novo, contribute to PUV formation, a male-limited phenotype.

TBK1 and TNFRSF13B mutations and an autoinflammatory disease in a child with lethal COVID-19.

Among children, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections are typically mild. Here, we describe the case of a 3.5-year-old girl with an unusually severe presentation of coronavirus disease (COVID-19). The child had an autoinflammatory disorder of unknown etiology, which had been treated using prednisolone and methotrexate, and her parents were half cousins of Turkish descent. After 5 days of nonspecific viral infection symptoms, tonic-clonic seizures occurred followed by acute cardiac insufficiency, multi-organ insufficiency, and ultimate death. Trio exome sequencing identified a homozygous splice-variant in the gene TBK1, and a homozygous missense variant in the gene TNFRSF13B. Heterozygous deleterious variants in the TBK1 gene have been associated with severe COVID-19, and the variant in the TNFRSF13B gene has been associated with common variable immunodeficiency (CVID). We suggest that the identified variants, the autoinflammatory disorder and its treatment, or a combination of these factors probably predisposed to lethal COVID-19 in the present case.

Mapping of cis-acting expression quantitative trait loci in human scalp hair follicles.

BACKGROUND: The association of molecular phenotypes, such as gene transcript levels, with human common genetic variation can help to improve our understanding of interindividual variability of tissue-specific gene regulation and its implications for disease. METHODS: With the aim to capture the spectrum of biological processes affected by regulatory common genetic variants (minor allele frequency ≥ 1%) in healthy hair follicles (HFs) from scalp tissue, we performed a genome-wide mapping of cis-acting expression quantitative trait loci (eQTLs) in plucked HFs, and applied these eQTLs to help further explain genomic findings for hair-related traits. RESULTS: We report 374 high-confidence eQTLs found in occipital scalp tissue, whose associated genes (eGenes) showed enrichments for metabolic, mitotic and immune processes, as well as responses to steroid hormones. We were able to replicate 68 of these associations in a smaller, independent dataset, in either frontal and/or occipital scalp tissue. Furthermore, we found three genomic regions overlapping reported genetic loci for hair shape and hair color. We found evidence to confirm the contributions of PADI3 to human variation in hair traits and suggest a novel potential candidate gene within known loci for androgenetic alopecia. CONCLUSIONS: Our study shows that an array of basic cellular functions relevant for hair growth are genetically regulated within the HF, and can be applied to aid the interpretation of interindividual variability on hair traits, as well as genetic findings for common hair disorders.

Sustained Immunoparalysis in Endotoxin-Tolerized Monocytic Cells.

Sepsis is associated with a strong inflammatory reaction triggering a complex and prolonged immune response. Septic patients have been shown to develop sustained immunosuppression due to a reduced responsiveness of leukocytes to pathogens. Changes in cellular metabolism of leukocytes have been linked to this phenomenon and contribute to the ongoing immunological derangement. However, the underlying mechanisms of these phenomena are incompletely understood. In cell culture models, we mimicked LPS tolerance conditions to provide evidence that epigenetic modifications account for monocyte metabolic changes which cause immune paralysis in restimulated septic monocytes. In detail, we observed differential methylation of CpG sites related to metabolic activity in human PBMCs 18 h after septic challenge. The examination of changes in immune function and metabolic pathways was performed in LPS-tolerized monocytic THP-1 cells. Passaged THP-1 cells, inheriting initial LPS challenge, presented with dysregulation of cytokine expression and oxygen consumption for up to 7 days after the initial LPS treatment. Proinflammatory cytokine concentrations of TNFα and IL1ß were significantly suppressed following a second LPS challenge (p < 0.001) on day 7 after first LPS stimulation. However, the analysis of cellular metabolism did not reveal any noteworthy alterations between tolerant and nontolerant THP-1 monocytes. No quantitative differences in ATP and NADH synthesis or participating enzymes of energy metabolism occurred. Our data demonstrate that the function and epigenetic modifications of septic and tolerized monocytes can be examined in vitro with the help of our LPS model. Changes in CpG site methylation and monocyte function point to a correlation between epigenetic modification in metabolic pathways and reduced monocyte function under postseptic conditions.

Advanced paternal age as a risk factor for neurodevelopmental disorders: a translational study.

Advanced paternal age (APA) is a risk factor for several neurodevelopmental disorders, including autism and schizophrenia. The potential mechanisms conferring this risk are poorly understood. Here, we show that the personality traits schizotypy and neuroticism correlated with paternal age in healthy subjects (N = 677). Paternal age was further positively associated with gray matter volume (VBM, N = 342) in the right prefrontal and the right medial temporal cortex. The integrity of fiber tracts (DTI, N = 222) connecting these two areas correlated positively with paternal age. Genome-wide methylation analysis in humans showed differential methylation in APA individuals, linking APA to epigenetic mechanisms. A corresponding phenotype was obtained in our rat model. APA rats displayed social-communication deficits and emitted fewer pro-social ultrasonic vocalizations compared to controls. They further showed repetitive and stereotyped patterns of behavior, together with higher anxiety during early development. At the neurobiological level, microRNAs miR-132 and miR-134 were both differentially regulated in rats and humans depending on APA. This study demonstrates associations between APA and social behaviors across species. They might be driven by changes in the expression of microRNAs and/or epigenetic changes regulating neuronal plasticity, leading to brain morphological changes and fronto-hippocampal connectivity, a network which has been implicated in social interaction.

Mutations in SREBF1, Encoding Sterol Regulatory Element Binding Transcription Factor 1, Cause Autosomal-Dominant IFAP Syndrome.

IFAP syndrome is a rare genetic disorder characterized by ichthyosis follicularis, atrichia, and photophobia. Previous research found that mutations in MBTPS2, encoding site-2-protease (S2P), underlie X-linked IFAP syndrome. The present report describes the identification via whole-exome sequencing of three heterozygous mutations in SREBF1 in 11 unrelated, ethnically diverse individuals with autosomal-dominant IFAP syndrome. SREBF1 encodes sterol regulatory element-binding protein 1 (SREBP1), which promotes the transcription of lipogenes involved in the biosynthesis of fatty acids and cholesterols. This process requires cleavage of SREBP1 by site-1-protease (S1P) and S2P and subsequent translocation into the nucleus where it binds to sterol regulatory elements (SRE). The three detected SREBF1 mutations caused substitution or deletion of residues 527, 528, and 530, which are crucial for S1P cleavage. In vitro investigation of SREBP1 variants demonstrated impaired S1P cleavage, which prohibited nuclear translocation of the transcriptionally active form of SREBP1. As a result, SREBP1 variants exhibited significantly lower transcriptional activity compared to the wild-type, as demonstrated via luciferase reporter assay. RNA sequencing of the scalp skin from IFAP-affected individuals revealed a dramatic reduction in transcript levels of low-density lipoprotein receptor (LDLR) and of keratin genes known to be expressed in the outer root sheath of hair follicles. An increased rate of in situ keratinocyte apoptosis, which might contribute to skin hyperkeratosis and hypotrichosis, was also detected in scalp samples from affected individuals. Together with previous research, the present findings suggest that SREBP signaling plays an essential role in epidermal differentiation, skin barrier formation, hair growth, and eye function.

A high-salt diet compromises antibacterial neutrophil responses through hormonal perturbation.

The Western diet is rich in salt, which poses various health risks. A high-salt diet (HSD) can stimulate immunity through the nuclear factor of activated T cells 5 (Nfat5)-signaling pathway, especially in the skin, where sodium is stored. The kidney medulla also accumulates sodium to build an osmotic gradient for water conservation. Here, we studied the effect of an HSD on the immune defense against uropathogenic E. coli-induced pyelonephritis, the most common kidney infection. Unexpectedly, pyelonephritis was aggravated in mice on an HSD by two mechanisms. First, on an HSD, sodium must be excreted; therefore, the kidney used urea instead to build the osmotic gradient. However, in contrast to sodium, urea suppressed the antibacterial functionality of neutrophils, the principal immune effectors against pyelonephritis. Second, the body excretes sodium by lowering mineralocorticoid production via suppressing aldosterone synthase. This caused an accumulation of aldosterone precursors with glucocorticoid functionality, which abolished the diurnal adrenocorticotropic hormone-driven glucocorticoid rhythm and compromised neutrophil development and antibacterial functionality systemically. Consistently, under an HSD, systemic Listeria monocytogenes infection was also aggravated in a glucocorticoid-dependent manner. Glucocorticoids directly induced Nfat5 expression, but pharmacological normalization of renal Nfat5 expression failed to restore the antibacterial defense. Last, healthy humans consuming an HSD for 1 week showed hyperglucocorticoidism and impaired antibacterial neutrophil function. In summary, an HSD suppresses intrarenal neutrophils Nfat5-independently by altering the local microenvironment and systemically by glucocorticoid-mediated immunosuppression. These findings argue against high-salt consumption during bacterial infections.

The role of environmental stress and DNA methylation in the longitudinal course of bipolar disorder.

BACKGROUND: Stressful life events influence the course of affective disorders, however, the mechanisms by which they bring about phenotypic change are not entirely known. METHODS: We explored the role of DNA methylation in response to recent stressful life events in a cohort of bipolar patients from the longitudinal PsyCourse study (n = 96). Peripheral blood DNA methylomes were profiled at two time points for over 850,000 methylation sites. The association between impact ratings of stressful life events and DNA methylation was assessed, first by interrogating methylation sites in the vicinity of candidate genes previously implicated in the stress response and, second, by conducting an exploratory epigenome-wide association analysis. Third, the association between epigenetic aging and change in stress and symptom measures over time was investigated. RESULTS: Investigation of methylation signatures over time revealed just over half of the CpG sites tested had an absolute difference in methylation of at least 1% over a 1-year period. Although not a single CpG site withstood correction for multiple testing, methylation at one site (cg15212455) was suggestively associated with stressful life events (p < 1.0 × 10-5). Epigenetic aging over a 1-year period was not associated with changes in stress or symptom measures. CONCLUSIONS: To the best of our knowledge, our study is the first to investigate epigenome-wide methylation across time in bipolar patients and in relation to recent, non-traumatic stressful life events. Limited and inconclusive evidence warrants future longitudinal investigations in larger samples of well-characterized bipolar patients to give a complete picture regarding the role of DNA methylation in the course of bipolar disorder.